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EXECUTIVE  SUMMARY 


We  measured  the  vapor-phase  absorptivity  coefficient  of  pinacolyl 
methylphosphonofluoridate  (GD,  Soman)  in  the  mid-infrared  (IR).  We  used  agent  filled 
saturator  cells  suspended  in  a  temperature  controlled  liquid  bath  to  generate  continuous  streams 
of  the  compound  diluted  in  nitrogen,  which  were  sent  to  a  variable  path  White  cell  and  measured 
using  a  high  resolution  research  grade  Fourier  transform  IR  spectrometer.  The  purity  of  the 
feedstock  was  verified  by  nuclear  magnetic  resonance  (NMR)  spectroscopy,  gas  chromatography 
(GC)-mass  spectrometry,  and  GC-IR.  The  mass  of  GD  in  the  vapor  was  determined  with  a 
gravimetric  method,  and  the  ratio  was  verified  with  thermal  desorption  GC.  Ten  spectra  at 
different  concentration-pathlength  products  were  processed  line  by  line  through  least  squares 
analysis  using  MatLab®  to  produce  the  absorptivity  coefficient  of  the  compound  and  the 
statistical  uncertainty  in  the  data.  Uncertainties  in  the  data,  expanded  to  a  confidence  interval  of 
2a  (P  =  0.95),  are  Type-A:  2.7%  and  Type-B:  2.8%  of  the  absorptivity  coefficient.  In  this  report, 
we  compare  our  results  to  the  data  obtained  by  another  laboratory  using  a  different  vapor 
generation  method. 
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PREFACE 


The  work  described  in  this  report  was  performed  under  the  direction  of  the 
Detection  Capability  Officer,  Defense  Threat  Reduction  Agency  Joint  Science  and  Technology 
Office.  This  work  was  started  and  completed  in  September  2008. 

The  use  of  either  trade  or  manufacturers’  names  in  this  report  does  not  constitute 
an  official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

This  report  has  been  approved  for  public  release.  Registered  users  should  request 
additional  copies  from  the  Defense  Technical  Information  Center;  unregistered  users  should 
direct  such  requests  to  the  National  Technical  Information  Service. 
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VAPOR-PHASE  INFRARED  ABSORPTIVITY  COEFFICIENT 
OF  PINACOLYL  METHYLPHOSPHONOFLUORIDATE 


1.  INTRODUCTION 

We  determined  the  high  resolution  vapor-phase  absorptivity  coefficients  of  the 
nerve  agent,  pinacolyl  methylphosphonofluoridate  (GD,  Soman)  in  the  spectral  range  of  4000- 
550  cm'1  in  units  of  (pmol/moiy'm  1  and  computed  the  uncertainties  in  the  data.  Previous  efforts 
at  the  U.S.  Army  Edgewood  Chemical  Biological  Center  (ECBC)  to  obtain  the  absorptivity 
coefficients  of  this  compound  date  back  to  the  1960s  and  were  performed  at  lower  resolution  on 
a  grating  spectrometer.  Given  the  extensive  use  of  vapor-phase  infrared  (IR)  spectra  for 
standoff  detection  of  chemical  warfare  agents  (CWAs),  the  need  for  current  data  is  apparent. 

Soman  has  the  Chemical  Abstracts  Service  (CAS)  Registry  Number  96-64-0  and 
is  indexed  under  the  name:  phosphonofluoridic  acid,  methyl-,  1,2,2-trimethylpropyl.  The 
molecular  formula  is  C7H16FO2P  and  it  has  a  molecular  weight  of  1 82. 1 78.  Other  synonyms  for 
the  compound  include  3,3-dimethyl-n-but-2-yl  methylphosphonofluoridate;  methyl  pinacolyl 
phosphonofluoridate;  pinacolyl  methylfluorophosphonate;  Zoman;  and  EA-I210.  The  structure 
is  shown  in  Figure  I. 


HUC- 


ch3  ch3 
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Figure  I .  Structure  of  Pinacolyl  Methylphosphonofluoridate  (GD) 

Pinacolyl  methylphosphonofluoridate  is  a  nerve  agent  that  was  first  synthesized  in 
Germany  in  1944  although  it  was  not  used  in  combat  during  WWII.  The  vapor  pressure  of  the 
compound  at  several  temperatures  is:  36.4  Pa  (20  °C),  53.4  Pa  (25  °C),  and  101325  Pa  (197.8  °C, 
calculated  boiling  temperature). 


2.  EXPERIMENTAL  PROCEDURE 

2.1  Instrumental  Details. 

The  system  used  to  generate  the  continuous  vapor  stream  was  an  adaptation  of  the 
saturator  cell  method  developed  at  ECBC  for  measuring  the  volatility  of  CWA  related 
compounds.3 


The  method,  modified  to  generate  continuous  streams  of  chemical  compounds  for 
obtaining  quantitative  vapor-phase  IR  spectra,  has  been  used  to  measure  the  absorptivity 
coefficients  of  benzene,4  as  well  as  a  variety  of  CWA  related  compounds.5  References  4  and  5 
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describe  the  experimental  setup  as  well  as  the  data  collection  and  post-processing  in  more  detail. 
The  saturator  passes  a  stream  of  nitrogen  carrier  gas,  obtained  from  the  boil-off  of  a  bulk  liquid 
nitrogen  tank,  across  a  cylindrical  alumina  wicking  mechanism  closed  at  one  end  in  a  glass 
holder  filled  with  the  analyte.  This  creates  a  saturated  vapor-liquid  equilibrium  of  the  analyte  on 
the  downstream  side  of  the  saturator  cell.  The  concentration  of  the  analyte  is  determined  by  the 
temperature  of  the  liquid  phase.  By  suspending  the  saturator  cell  in  a  constant  temperature  bath, 
the  concentration  of  the  analyte  can  be  predicted  by  its  vapor  pressure  at  the  temperature  of  the 
bath.  The  Quantitative  FTIR  Laboratory  is  equipped  with  a  Brooks  (Brooks  Instrument  Co., 
Hatfield,  PA)  Model  5850S  mass  flow  controller  to  maintain  a  constant  flow  to  the  saturator  cell 
along  with  a  second  mass  flow  controller  to  add  diluent  to  the  stream,  providing  an  additional 
means  of  adjusting  the  concentration  of  the  compound  delivered  to  the  White  cell  of  the  FTIR. 
Linearity  of  the  S  series  mass  flow  controllers  is  adjusted  using  a  second  order  polynomial, 
resulting  in  uncertainties  of  approximately  1%  or  better  of  the  flow  rate  at  rates  >25%  of  full 
scale. 


Spectra  were  obtained  with  a  Bruker  (Bruker  Optics,  Billerica,  MA)  Model 
1FS/66V  Fourier  transform  1R  spectrometer  (FTIR).  The  instrument  is  equipped  with  both 
deuterated  triglycine  sulfide  and  mercury-cadmium-telluride  (HgCdTe)  detectors  and  is  capable 
of  obtaining  spectra  with  a  maximum  spectral  resolution  of  0. 1 1 25  cm  1  (unapodized).  The 
interferograms  were  recorded  from  15798-0  cm  1  with  a  resolution  of  0.125  cm  Absorbance 
(log  base-1 0)  spectra  were  processed  with  boxcar  apodization  and  2X  zero  filled  to  obtain  a  data 
spacing  of  0.0625  cm"1.  The  instrument  is  equipped  with  a  variable  path  White  cell.  The 
experimental  data  used  path  lengths  of  2.727  and  5.377-m.  The  temperature  of  the  White  cell 
was  maintained  at  23  ±  0.1  °C  by  using  a  thermostatically  controlled  chamber  enclosing  the 
spectrometer  and  cell.  Data  were  acquired  at  a  speed  of  60  KHz  (HeNe  laser  zero  crossing 
frequency)  using  the  HgCdTe  detector.  Single  beam  spectra  of  the  CWA  were  ratioed  against 
spectra  of  clean,  dry  nitrogen.  To  minimize  the  effects  of  nonlinearity  in  the  detector,  the 
interferograms  were  processed  using  the  proprietary  Opus  nonlinearity  correction  function.  All 
interferograms  were  archived  enabling  further  post-processing  of  data. 

Temperature  and  pressure  data  were  recorded  using  National  Institute  of 
Standards  and  Technology  traceable  digital  manometers  and  thermometers,  and  all  the  data  were 
archived.  Concentration-pathlength  products  (CL)  were  computed  in  units  of  pmol/mol(m) 
(ppm-m).  A  differential  pressure  manometer  had  previously  been  used  to  measure  the  dynamic 
pressure  in  the  White  cell  with  gas  flowing  into  the  cell,  and  the  ambient  pressure  was  plotted 
versus  the  differential  pressure.  The  resulting  equation  was  used  to  correct  the  readings  from  the 
ambient  pressure  manometer  to  the  pressure  in  the  White  cell.  The  concentration-pathlength  data 
were  corrected  to  296  K  and  1.0132  X  10'  Pa  (760  Torr)  using  the  ideal  gas  law. 
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2.2 


Feedstock. 


The  material  used  to  generate  the  vapor  streams  for  the  experiments  was  purified 
by  fractional  distillation  to  remove  the  stabilizer.  The  purity  was  determined  to  be 
1 0 1 .6  ±  1.4  mol%  3 1 P  nuclear  magnetic  resonance  spectroscopy  (NMR)  using  an  internal 
standard;  and  99.8%  gas  chromatography-mass  spectrometry  (GC-MS),  with  0.2%  pinacolyl 
alcohol.  Analysis  by  GC-1R  confirmed  the  presence  of  the  pinacolyl  alcohol  in  the  feedstock. 
Sampling  of  the  vapor  effluent  by  thermal  desorption-GC  with  a  flame  ionization  detector 
indicated  a  purity  of  99.0%.  The  spectra  of  the  vapor  from  the  saturator  cell  showed  the 
presence  of  both  water  vapor  and  hydrogen  fluoride  in  the  effluent  at  a  combined  mass 
concentration  that  ranged  from  0.25  to  0.5%.  We  used  a  concentration  of  99.5  ±  1.0  (2a)  mass% 
to  compute  the  concentration-path  length  products  of  the  vapor. 


3.  RESULTS  AND  DISCUSSION 

Three  trials  were  run  to  obtain  spectra  at  1 1  concentration-pathlength  products.  A 
trial  is  defined  as  filling  and  weighing  the  saturator  cell,  suspending  it  in  the  bath,  applying  a 
stream  of  nitrogen  for  a  measured  time,  acquiring  several  spectra,  stopping  the  nitrogen  and 
removing  it  from  the  bath,  and  reweighing  the  saturator  cell  after  drying  the  exterior  surfaces  and 
re-equilibrating  to  room  temperature.  The  constant  temperature  bath  was  maintained  at  23.0  °C 
for  all  three  experiments. 

Trace  water  vapor  and  hydrogen  chloride  were  observed  in  several  spectra.  The 
concentrations  of  the  compounds  in  the  vapor,  which  decreased  over  time,  were  computed  using 
their  absorptivity  coefficients  and  the  features  of  the  impurities  were  removed  by  spectral 
subtraction.  Baseline  corrections  of  the  spectra  were  performed  with  a  two  point  linear 
subtraction  using  a  MatLab  algorithm  written  in-house.  The  baseline  correction  in  no  case 
exceeded  0.001  A[A  =  -logio(T)]. 

The  composite  spectrum  (absorptivity  coefficient)  was  computed  using  spectra 
with  concentration-path  length  products  ranging  50-416  pmol/mol(m)  (corrected  to  296  K  and 
101325  Pa).  As  an  initial  check  of  the  quality  of  the  data,  Beer’s  Law  plots  of  two  spectral  lines 
(1212.5  and  720.6  cm'1)  were  calculated  using  MatLab.  These  plots  showed  that,  at  least  for 
these  two  spectral  lines,  the  data  appeared  to  be  well  fitted  with  no  points  lying  outside  the  95% 
confidence  limits  for  a  repeated  set,  a  repeated  single  x,  or  the  95%  confidence  limits  for  a 
Grubbs  Test  for  Outliers6  (Figures  2  and  3). 
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Figure  2.  Beer’s  Law  Plot  of  1325.97  cm  Line  in  Vapor-Phase  Spectrum  of  GD,  Soman 
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Figure  3.  Beer’s  Law  Plot  of  1020.49  cm  1  Line  in  Vapor-Phase  Spectrum  of  GD,  Soman 
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The  absorptivity  coefficient  (a)  and  uncertainty  (Type-A,  2a)  were  computed  line 
by  line  within  the  spectral  range  of  4000-550  cm'1  using  a  MatLab  program  written  in-house. 
Values  of  ( A  =  -log7)  >  1 .5  are  normally  assigned  a  weight  of  zero.  Because  all  values  recorded 
were  for  A  <  1 .368,  all  data  were  weighted  at  1 .  Figure  4  is  the  plotted  absorptivity  coefficient 
(a),  and  Type-A  uncertainty  for  the  computed  spectral  range.  The  figures  are  plotted  with  a  in 
(pmol/moiy'm1.  To  obtain  a  in  (mg/m2)'1,  multiply  the  values  in  the  mantissa  of  Figure  4  by 
0. 1 326.  This  factor  is  derived  from  eq  1  using  the  molecular  weight  of  Lewisite  (182.1 78). 

m2  (  24,15 

mg  ^  mw 


\  _  mol 
)  /.imol(m) 


Figure  4.  Absorptivity  Coefficient  and  Statistical  Uncertainty  (Type-A,  2a)  for  GD,  Soman 

Table  1  provides  the  absorptivity  coefficients  in  (pmol/mol)  'm  1  and  (mg/m2)'1 
for  selected  bands  in  units  of  wavenumber  and  micrometers  (pm). 


Table  1.  Absorptivity  Coefficient  of  Pinacolyl  Methylphosphonofluoridate  for  Selected  Bands 


P7(cm  ')  2975.29 

A/pm  3.361 

(pmol/moiy’m  1  8.314  X  1  O'04 

m2/mg  1.102  XI  O'04 


1325.95 

7.542 

1 .054  X  10'03 
1.397  X  10'04 


1020.47 

9.799 

3.274  X  10'03 
4.340  X  1  O'04 


921.69 
10.850 
7.205  X  10'04 
9.551  X  10'05 
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In  general,  expanded  Type-A  uncertainties  were  1-1.5%  of  the  absorptivity 
coefficient  as  seen  in  Figures  5  and  6.  Figure  5  is  a  plot  of  absorptivity  coefficients  (abscissa) 
and  fractional  uncertainty  (Type-A,  Ua,  2o)  (mantissa). 


Figure  5.  Absorption  Coefficient  (abscissa)  and  Type-A  Uncertainty  (Fractional,  2o)  for  Soman 


Figure  6  is  a  plot  of  the  absorptivity  coefficient  and  uncertainty.  It  includes  a  best 
fit  of  the  data  points  obtained  by  least  squares,  which  is  an  approximation  of  Ua-  ~  ax+b.  For  the 
fitted  line  in  Figure  6,  the  coefficients  are  ,4  =  1.24  X  l(f2  and  b=  1.95  X  10'6  (Table  2). 
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Figure  6.  Absorption  Coefficient  (abscissa)  and  Type  A  Uncertainty  (2a)  for  GD,  Soman. 
The  line  indicates  the  “best  fit”  obtained  by  least  squares. 


Table  2.  Type-A  Statistical  Uncertainty  for  GD.  Soman  Vapor-Phase  Absorptivity  Coefficient 


Type-A  Uncertainty 

2a  ~  ma  +  h 

Slope 

Intercept 

m 

b 

1.24  X  10‘2 

1.95  X  10'6 

Type-B  estimated  standard  errors,  along  with  their  sources,  as  well  as  the 
combined  Type-A  and  B  uncertainties  are  provided  in  Table  3.  The  expanded  combined  Type-B 
uncertainty  was  computed  using  eq  2: 

Ab  =  ( A  L2+ AT2+ A  P2+ A  FT  I R2+ AN  L2+A  M  R2) 1  /2  (2) 


The  sources  of  uncertainty  and  their  fractional  values,  as  well  as  an  explanation  of 
the  symbols  in  eq  1  are  given  in  Table  4.  The  Type-A  uncertainty  of  2.7%  (2a)  indicated  in 
Table  3  is  valid  for  values  of  the  absorptivity  coefficient  >0.000133  (pmol/mol)''m  1  (a>5%  of 
the  peak  of  maximum  intensity).  For  the  most  intense  features,  the  Type-A  uncertainty  averages 
1 .2%  of  the  absorptivity  coefficient.  Among  the  Type  B  uncertainties,  the  detector  nonlinearity 
dominates  at  0.01  (la). 
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Table  3.  Uncertainties  in  Absorptivity  Coefficient  of  GD,  Soman  from  ECBC  Data  where 
a  >  0.000164  (pmol/moiy'm"1 


Symbol 

Fractional  Deviation 

Source 

AL 

0.005 

Pathlength 

AT 

0.0006 

Temperature  of  White  cell 

AP 

0.002 

Pressure 

AFTIR 

0.0005 

Drift  in  spectrometer 

ANL 

0.01 

Nonlinearity  in  detector 

AMR 

0.005 

Mass  rate 

AD 

0.005 

Dilution  rate 

Apurity 

0.005 

Purity  of  feedstock 

Ar 

0.028 

Combined  type  B  (2a) 

Aa 

0.027 

Type-A  deviation  (2a) 

4.  INTERLABORATORY  COMPARISON 

Comparison  of  data  between  laboratories,  especially  when  obtained  using 
different  methods,  can  be  useful  in  assessing  the  experimental  methods  and  for  validating  the 
accuracy  of  the  data.  One  other  source  for  comparing  the  high  resolution  absorptivity  coefficient 
of  pinacolyl  methylphosphonofluoridate  to  the  results  acquired  at  ECBC  are  the  data  obtained  by 
Pacific  Northwest  National  Laboratory  (PNNL)  at  Dugway  Proving  Ground.  The  PNNL 
database  includes  a  detailed  discussion  of  experimental  methods  and  statistical  processing.  Each 
absorptivity  coefficient  spectrum  is  provided  with  a  metadata  file  listing  physical  data  for  the 
compound,  temperature  and  pressure  range,  number  of  individual  spectra  included  in  the  fit  of 
the  absorptivity  coefficient,  and  Type-A  and  Type-B  uncertainties. 

As  can  be  seen  in  Figure  7,  the  spectra  from  ECBC  and  PNNL  are  qualitatively 
similar.  The  metadata  file  associated  with  the  PNNL  spectrum  indicates  uncertainties  of  2.0% 
(Type-A)  and  <10%  (Type-B).  A  comparison  of  the  integrated  absorptivity  coefficients  from  the 
two  laboratories  shown  in  Table  4,  however,  shows  that  the  differences  between  the  spectra 
appear  to  be  larger  than  their  uncertainties.  As  seen  in  Figure  8,  the  baselines  of  the  two  spectra 
differ  significantly.  In  general,  across  the  spectral  range,  the  baseline  of  the  PNNL  spectrum 
appears  to  be  shifted  on  the  lower  frequency  (“red”)  side  of  absorption  features.  While  this  may 
in  part  arise  from  the  third  order  polynomial  baseline  subtraction  used  to  correct  the  PNNL 
spectrum,  we  believe  this  is  more  likely  the  effect  of  scattering  from  aerosols  present  in  the 
vapor  of  the  spectra  to  compute  the  composite  spectrum.  The  strong  band  associated  with  P  =  O 
stretch  appears  in  GD  at  1301  cm  1  in  the  vapor-phase  and  1287  cm'1  in  the  liquid.  Thus  the 
simplest  method  to  detect  the  presence  of  liquid  in  the  vapor  is  to  subtract  an  aerosol-free 
spectrum  from  the  PNNL  spectrum  and  compare  the  residuals  in  that  region.  Indeed,  when  the 
ECBC  spectrum  is  subtracted  from  the  PNNL  spectrum,  the  residual  (Figure  9)  does  show  the 
presence  of  absorption  features  apparently  associated  with  liquid-phase  material  in  the  vicinity  of 
1300  cm'1.  The  aerosols  are  likely  caused  by  the  syringe  pump  method  used  to  generate  the 
vapor  used  for  the  PNNL  spectra.  The  apparent  presence  of  the  liquid  in  the  PNNL  spectrum. 
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combined  with  the  effects  of  aerosol  scattering  effects  on  the  baseline,  renders  a  quantitative 
comparison  of  the  vapor  spectra  from  the  two  laboratories  difficult  at  best. 


x  10'3 


Figure  7.  Absorptivity  Coefficient  of  Vapor-Phase  Soman  from  PNNL.  The  spectrum  is 
qualitatively  similar  to  the  spectrum  from  our  laboratory  as  seen  in  Figure  4. 


Table  4.  Comparison  of  Integrated  Absorptivity  Coefficients  of  GD,  Soman  from  PNNL  and 
ECBC 


Integrated  Absorptivity 
Coefficient 

Range/cm  1 

PNNL 

ECBC 

(ECBC-PNNL)/ECBC  X  100 

4000-550 

0.2791 

0.2476 

-12.7% 

1540-683 

0.2131 

0.1990 

-7.1% 

3258-2423 

0.0593 

0.0459 

-29.2% 
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X  id'4 


Wavenumber  (cm'1) 


Figure  8.  Absorptivity  Coefficient  Spectra  in  Vicinity  of  3000  cm  1  from  ECBC  and  PNNL.  The 
PNNL  spectrum  appears  to  exhibit  a  baseline  shift  on  the  lower  frequency  side  of  the  peaks 
associated  with  CH  stretch.  This  may  be  indicative  of  aerosol  scattering  in  the  PNNL  spectrum. 
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x  10'4 


Wavenumber  (cm'1) 

Figure  9.  Absorptivity  Coefficient  Spectra  of  GD  from  ECBC  (vapor)  and  PNNL  (liquid)  and 
Residual  Spectum  Obtained  by  Subtracting  PNNL  Vapor  Spectrum  from  ECBC  Vapor 
Spectrum.  Arrows  indicate  the  presence  of  absorption  features  in  the  vicinity  of 
1300  cm  1  that  are  apparently  associated  with  liquid-phase  GD  present  in  the  vapor  used  to 
generate  the  PNNL  vapor-phase  spectrum. 


5.  CONCLUSIONS 

We  used  a  saturator  cell  system  to  acquire  vapor-phase  spectra  of  the  nerve  agent, 
pinacolyl  methylphosphonofluoridate.  These  spectra  were  then  used  to  compute  the  vapor-phase 
absorptivity  coefficient  at  0.125  cm'1  resolution.  The  uncertainties  are  Type-A  =  2.7%  and 
Type-B  =  2.8%  of  the  absorptivity  coefficient  for  absorptivity  coefficient  >5%  of  the  peak  of 
maximum  intensity. 
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